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Abstract. Recent years have been marked by substantial changes in our understanding 
of magnetic turbulence and magnetic reconnection, which, in its turn induced better under- 
standing of cosmic ray diffusion and acceleration. Current models of magnetized turbulence 
are no more ad hoc constructions, but numerically tested theories. In this very short review 
we summarize topics presented in two talks given at the conference and provide a brief 
sketch of the vast and rapidly developing field. We discuss how turbulence decreases the 
efficient mean free path of the particles in the collisionless plasmas in galaxy clusters and 
claim that this makes MHD turbulence description applicable to a wider range of scales. We 
discuss the properties of MHD turbulence and its relation to magnetic reconnection. Finally, 
we overview how turbulence induces particle acceleration via second order Fermi process 
and affects first order Fermi acceleration in shocks and reconnection regions. 
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1 . Guide to the review 

Mergers between galaxy clusters are the 
most energetic events in the present day 
Universe. During these mergers a fraction 
of the gravitational energy can be con- 
verted into fluid motions, i.e. shocks and 
turbulence, that generate magnetic fields 
and, through a variety of processes, ac- 
celerate ^dativistic_OTotons__and_decfiom 
(e.g., iRvu et alTl l2003t ICassano & Brunettil 



20051; iBrunetti & Lazarian I 120071 
Hoeft & Brueg genl 120071; Pfromme r et al. 
2008; Skillma neTal. I bOOSMBrunetti et al" 
20091: IVazza et al. 1 120091) . In this short review 



we address some of the basic processes in- 



volved, namely, magnetic turbulence in galaxy 
clusters and the possibility of its observational 
studies (Sect. 2), magnetic reconnection (Sect. 
3), as well as various ways of accelerating cos- 
mic rays (Sect. 4). Our summary is presented 
in Sect. 5. 

2. Turbulence in clusters of galaxies 

Astrophysical fluids are characterized by 
high Reynolds numbers and are known 
to be turb ulent (e.g., [A rmstron g et alH 



[19951 

CI m cti~r 



20 lOt 
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Chepurnov & Lazarian 
Elme green & Scalo H2004tlMcKee & Ostrikerl 
2007). As properties of turbulent magnetized 
fluids are very different from laminar ones, the 
correct description of the particle acceleration 
requires taking into account the fundamental 
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properties of magnetic turbulence as well as 
the mutual feedback of magnetic fields and 
cosmic rays in the turbulent fluids. 

2.1. Properties of intracluster plasmas: 
instabilities induced by turbulence 
and effective collisions 

Turbulence in galaxy clusters is magnetized. 
A very important question is whether the 
MHD description of turbulence is applica- 
ble. When Coulomb collisions in the rarefied 
inter-galactic medium (IGM) are considered 
one has to conclude that the plasma is col- 
lisionless. This strongly affects the prepara- 
tion of compressibl e modes, cosmic ray ac- 
celeration etc (see lBrunetti & Lazarian l d20Q7h 
and ref. therein). In what follows we ar- 
gue that the degree of collisionality of astro- 
physical plasmas is underestimated when only 
Coloumb collis i ons are ta ken into account (see 
lLazarian et al. I d2010l) : iBrunetti & Lazarianl 
1 201 lah \ 

It is well known that the mean free path 
of thermal protons due to Coulomb collisions 
in the h ot IGM is very large, ten to hundred 
kpc (e.g. JSarazin I1986I) . Fluids in such a colli- 
sionless regime can be ver y different from their 
collis i onal c ounte rparts dSchekochihin et alTl 
I2005L 120061 12010b . Several instabilities (e.g. 
firehose, mirror, gyroresonance etc) can be 
generated in the IGM in the presence of tur- 
bulence, leading to a transfer of the energy of 
large-scale compressions to perturbations on 
smaller scales. 

Many instabilities have growth rate which 
peaks at scales near the particle gyroradius, 
making very large the scale separation between 
the energy injection scale and the scale where 
this energy is being deposited. The scattering 
induced by instabilities dramatically reduces 
the effective mean free path of thermal ions 
decreasing the effective viscosity of the IGM 
and making plasmas effectively collisional on 
smaller scales. Indeed, charged particles can 
be randomized if they interact with perturbed 
magnetic field. If this field is a result of plasma 
instabilities, the process can be viewed as the 
collective interaction of an individual ion with 



the rest of the plasma, which is the process 
mediated by magnetic field. As a result, the 
fluid would behave as collisional on scales less 
that the Coulomb mea n free path. This issue 
has be en addressed in Lazarian & Beresnvak 
(2006) for the case of a collisionless fluid sub- 
ject to the gyroresonance instability that is 
driven by the anisotropy of the particle dis- 
tribution in the momentum space that arises 
from magnetic field compression; the larger 
the magnetic field compression, the higher the 
anisotropy induced and the higher is the insta- 
bility growth rate. They found that the turbu- 
lent magnetic compressions on the scale of the 
mean free path and less are the most effective 
for inducing the instability^. As the scattering 
happens on magnetic perturbations induced by 
the instability, the mean free path of particles 
decreases as a result of the operation of the in- 
stability. This results in the process being self- 
regulating, i.e. the stronger the turbulence at 
the scale of injection, the smaller is the mean 
free path of plasma particles and the larger is 
the span of scales over which the fluid behaves 
as essentially collisional. 

This induces an interesting picture where 
the mean free path of plasma protons depends 
on the level of compressions induced by tur- 
bulence and the mean free path is determined 
not by Columb collisions, but scattering on 
magnetic field inhomogeneities at the Larmor 
radius of thermal protons. The peculiar fea- 
ture of this picture is that the aforementioned 
magnetic field perturbations are not part of 
the normal turbulent cascade, but results of 
compressible turbulent motions at much larger 
scales. Thermal protons do not scatter each 
other through electric interactions, but partic- 
ipate in non-local interactions mediated by the 
perturbed magnetic field. The higher the level 
of comressible turbulence, the better is MHD 
description of the IGM. 



1 The larger scale compressions do still induce the 
instability, but their effect is reduced due to their re- 
duced ability to induce large changes of B over the 
time scale between scatterin g. The model is furthe r 
elaborated and improved in Yan & Lazarianl d201l|) . 
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2.2. MHD turbulence: brief summary of 
theory and main properties of 
turbulence in the IGM 

The last decade has been marked by 
substantial advances in understanding 

the MHD 



regime (e.g., Goldreich & Sridhar 


1995; Lazarian & Vishniac 1 


1999; 


Cho & Vish 


niacl 2000: Muller&Bis 


camrj 


2000; 


Lithwick & Goldreich 


200 1[ 


Cho et all 120021 Cho & Lazarian 1 


20021 


2003; Beresnvak & Lazarian 


2010; 


Kowal & Lazarian 2010). 



The presence 
makes MHD 
(Montgomery & Turner 



of a magnetic field 
turbulence aniso tropic 



S3; 



1981 



1984 



1983b Iffigdonl 
The relative importance 



Matthaeus et al. 
Oughtonl 120031) . 
of hydrodynamic and magnetic forces changes 
with scale, so the anisotropy of MHD tur- 
bulence does too. A lan dmark event in this 
was a seminal work by Goldreich & Sridharl 
d 1995b (GS95) which contained ideas that 
radically changed the further development of 
the subject. The corner stone of this model 
was the so-called critical balance idea which 
provided the analytical relation between the 
fluctuations parallel and perpendicular to 
the magnetic field. It also contains prophetic 
statements about mode coupling, providing 
guidelines for generalization of the model 
from the incompressible to compressible 
MHD. 

The original model was improved in the 
subsequent publications. For instance, GS95 
uses the closure relations that employ in the 
global system of reference related to the mean 
field, which, in fact, is an incorrect sys- 
tem to be used for the critical balance de- 
scription. In lLazarian & Vishniac 1 d!999h and 
later publications (e.g ., ICho& Vishniac f 20001 ; 
iMaron & Goldreich 1 1200 ll ICho et al. I 120021) 
the importance of the local system of reference, 
which is defined by the local direction of the 
magnetic field of a wave packet, was reviled. 
The local system of reference was employed 
in the su ccessful testing of the GS95 model. In 
addition. lLazarian & Vish niac l (fl999l) general- 
ized the GS95 model for the case when the tur- 



bulent injection velocity at the injection scale 
is less than the Alfvenic velocity. 

The predictions of the GS95 model 
are in rough agr eement with nu merical 
|Ch o & Vishniac 



simulations (e.g.. 



I200incho et al 



2000; 
2002; 



Maron & Goldreich 
Beresnvak & Lazaria n 2006), although some 
disagreement in terms of the measured spectral 
slope was noted. This disagreement produced 
a flow of papers with suggestions to improve 
the GS95 model by including additional 
effect s like d ynamical alignment dBoldvrev I 
I2005L 20061). polari z ation intermittency 
(Beresnvak & Lazarian 1 20061) . non-locality 
dGogo beridze et al. 2007 j). Mo r e rece n t stud- 
les in lBeresnvak& Lazarian I d2009L [2010) 
indicate that numerical simulations do not 
have sufficiently extended inertial range to get 
the actual spectral slop^l and therefore worries 
about the "inconsistency" of the GS95 model 
are premature. Evidence of the GS95 spectrum 
for the MHD inco mpressible turbulen ce was 
recently obtained bv lBeresnvakl (1201 11) . 

We shall add parenthetically that in a 
number of applications the empirical so-called 
composite 2D/slab model of magnetic fluctua- 
tions is used. In the latter model, which is also 
known as two-component model, it is assumed 
that fluctuations can be described as a superpo- 
sition of fluctuations with wave vectors parallel 
to the ambient large-scale magnetic field (so- 
called slab modes) and perpendicular to the 
mean field (so-called two-dimensional modes). 
It results in a maltese cross structure of mag- 
netic correlations. This model was developed 
to account for the solar wind observations, 
which it does well by adjusti ng the intensity 
of the two components (e.g. iMatthaeus et al. I 
1990). This theory of 2D fluctuations is 
consistent with th e theory of weak Alfvenic 
turbulence (e.g., Ng & Bhattachariee 19961; 



Lazarian & Vis hniac I 119991 Galtie r et al. I 
20001) but it can describe Alfvenic turbu- 



lence only over a limited range of scales. 



2 iBeres nvak & L azarian noticed that the 

magnetic turbulence is less local compared with the 
hydrodynamic one and therefore one requires a sub- 
stantially larger resolution to distinguish the actual 
spectral slope from the slope affected by the bottle- 
neck. 
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It may be treated as a parameterization of 
a particular type of magnetic perturbation 
dominated by the peculiarities of driving, but 
recent simulations by Gosh (2011, private 
communication) show that, at best, the model 
represents a special transient state of a not 
fully developed turbulence. In addition, slab 
modes do not arise naturally in turbulence 
with large-scale drivi ng, as shown by MHD 
nume rical simulations dCho & Lazarian Il2002l 
2003). Thus we do not consider this model for 
clusters of galaxies. 

The GS95 model of turbulence can be 
adopted to describe the Alfvenic part of 
MHD turbulent fluctuations in galaxy clus- 
ters. The model can be generalized also 
to compressible turbulence and even for 
supersonic motions numerical calculations 
show that the Alfvenic perturbations exhibit 
GS95 scaling (ICho & Lazarian 1 120021 l20()l 
IKowal & Lazarian 2010h . Wenote that we 
consider MHD turbulence where the flows of 
energy in the opposite directions are balanced. 
When this is not true, i.e. when the turbulence 
has non-zero cross-helicity, the properties of 
turbulence depart substantially from the GS95 
modefl. Solar wind presents a system with im- 
balanced turbulence. However, the degree of 
imbalance of turbulence in clusters of galaxies 
is unclear and we know that in compressible 
media the imbalance decreases due to reflect- 
ing of waves from pre-existing density fluc- 
tuations and due to t he development of para - 
metric instabililites Zann a et al. 2001). 
Similarly, we shall not discuss MHD turbu- 
lence at high magnetic Prandtl numbers, when 
the viscosity is much larg er than resistivity 
(e.g. JCho et aO2002Ll2003h . 

The GS95 model of turbulence combined 
with several considerations on the macro- and 
micro-physics of the IGM allows for a basic 
picture of the properties of turbulence in g alaxy 
clusters (e.g., Bru netti & Lazarian 2003). 



3 Among the existing theories of imbal- 
anced turbulence (e .g., iLithwick et al. I 120071: 

IChandranl 
all, but 



Beresnvak & Lazarian 
20081; iPerez & Boldyrev] 12009 



Beresnvak & Lazarian 




Cosmological numerical simulations show 
that large-scale turbulent motions are gener- 
ated during t he pro ce ss of cluster formation 



(Dolag et al 



. (120051) : llapichino & Niemeverl 
(2008): IVazzaetal.l (1201 lh . see also Nagai 
201 1, Iapichino 201 1, Vazza 201 1, this confer- 
ence). These motions, injected at large scales 
L„ ~ 300 - 500 kpc, are believed to provide 
the driver for turbulence at smaller scales. The 
typical velocity of the turbulent eddies at the 
injection scale is expected to be around Vl ~ 
500 - 700 km/s which makes turbulence sub- 
sonic, but strongly super- Alfvenic. Turbulence 
at large scales is thus essentially hydrodynamic 
and - most likely - made of a mix of compres- 
sive and incompressive eddies. The cascading 
of compressive (magnetosonic) modes may in- 
deed couple with that of solenoidal motions 
(Kolmogorov eddies). 

Viscosity in a turbulent and magnetised IGM 
is strongly suppressed due to the effect of the 
bending of magnetic field lines and of the 
perturbations of the magnetic field induced 
by plasma instabilities (e.g., Sect. 2.1). The 
important consequence is that an inertial range 
in the IGM is established - for both solenoidal 
and compressive modes - down to collisionless 
scales where a fraction of the turbulent energy 
is channelled into accele ration/heating of CR 
and th e rmal pl asma (see iBrunetti & Lazarian I 
d2007l 1201 lab for discussion). At small 
scales - in the inertial range - the velocity of 
turbulent eddies becomes sub-Alfvenic and 
turbulence is described in the MHD regime. 
At these scales the coupling between Alfven 
and compressible modes gets changed and 
only slow mod es are cascaded by A l fvenic 



modes (e.g.. IGoldreich & Sridhar I 119951: 



contradict to 



nume rical testing in [Beresnvak & Lazarian | d2009l 
2010) 



I Lithwick & Goldreich 2001; Cho & Lazarian I 
2002). The cascading of fast modes is not 
particularly sensitive to the presence of the 
other modes, fast modes remain isotropic 
while the spectrum of other modes becomes 
anisotropic. 

2.3. Spectroscopic ways of turbulence 
studies 

Recent observational advances to constrain tur- 
bulence in the IGM focussed on the broaden- 



Lazarian & Brunetti: Turbulence Reconnection and CRs 



7 



ing of lines in the X-ray spectra of galaxy clus- 
ters and provide interes ting limits in the case of 
cool-core clusters (e.g.. San ders et al71l201 ll) . 

Turbulence in clusters of galaxies can be 
studied in future using Doppler broderned 
emission. Here we briefly review techniques 
originally developed for studies of Doppler 
broderned emission and absorption lines in 
the interstellar medium research. These tech- 
niques, Velocity Channel Analysis (VCA) 
and Velocity Correlati on Spectrum (VCS) 
have been developed by Lazarian & Pogosvanl 
(2000, 12004 H006, 2008)) (henceforth LP00, 
LP04, LP06, LP08, respectively) and suc- 
cessfully used for studying turbulence in 
diffuse and molecu l ar gas (iLazarian l 2009t 



iPadoanet all 120091; IChepurnov et"aTTl2010b 



These techniques can be applied - at some ex- 
tent - to the case of the IGM and future X-ray 
telescopes with very good spectroscopic ca- 
pabilities (eg ASTRO-H) can be used for the 
studies. 

The difference between the VCA and the 
VCS is how the data is being handled. 
In the VCA technique the Position-Position- 
Velocity data cubes available through spec- 
troscopic observations are analysed by taking 
spectrum of the velocity slice of the cube. The 
spectrum of the fluctuations is analysed while 
changing the thickness of the velocity slice and 
the analytical description of the statistics of 
the fluctuations in the PPV slices described in 
LP00 and LP04 is used to obtain both the spec- 
trum of velocity and the spectrum of density 
fluctuations. 

A different approach is used in the VCS tech- 
nique, where fluctuations are analysed along 
the velocity coordinate. For the VCS technique 
one does not require good coverage of the 
Position-Position plane and a few spectral lines 
are sufficient to get the spectra of velocity and 
density (see Figure[TJ. 

New effects arise when strong absorption 
lines, which are in a saturated regime, are stud- 
ied. The procedure for studying of the saturated 
lines is presented in LP08. 

Our study of the effect of finite temper- 
atures for the technique reveals that, unlike 
the VCA, the temperature broadening does not 
prevent the turbulence spectrum from being re- 



Cluster 
of stars 




Effective beam 



Fig. 1. Illustration of VCS absorption studies of 
turbulence. Upper Panel: Schematic of measur- 
ing turbulence with absorption lines from point 
sources, e.g. stars, and an extended source, e.g. a 
galaxy. Lower panel: Velocity Coordinate Spectrum 
obtained using sampling of a turbulent volume 
along 10 lines of sight. The solid line corresponds 
to the theoretical expectatio ns. Readapted from 
Chep urnov & Lazarian 1l l2009h . 



covered from observations. Indeed, in VCA, 
gas temperature acts in the same way as the 
width of a channel. Within the VCS the term 
with temperature gets factorized and it influ- 
ences the amplitude of fluctuations (LP06). 
One can correct for this terrrQ, which also al- 
lows for a new way of estimating the interstel- 
lar gas temperature. 

Another advantage of the VCS compared to the 
VCA is that it reveals the spectrum of turbu- 
lence directly, while within the VCA the slope 
of the spectrum should be inferred from vary- 



4 To do this, one may attempt to fit for the tem- 
perature that would remove the exponential fall off 
in the spe ctrum of fluctuations along the velocity co- 
ordinate ( Chepu rnov & Lazarian 2006) 
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ing the thickness of the channel. As the thermal 
line width acts in a similar way as the channel 
thickness, additional care (see LP04) should be 
exercised not to confuse the channel that is still 
thick due to thermal velocity broadening with 
the channel that shows the thin slice asymp- 
totics. A simultaneous use of the VCA and the 
VCS makes the turbulence spectrum identifica- 
tion more reliable. 

Both VCA and VCS are applicable to stud- 
ies of not only emission, but also absorption 
lines. We note, that while dealing with emis- 
sion lines we may face additional complica- 
tions. For instance, Lazarian & Pogosyan (see 
LPOO, LP04, LP06, LP08) treated the emis- 
sivities proportional to the density to the first 
power. Therefore, in terms of scalings, the 
emissivities and densities were interchange- 
able. This is not true, however, when the emis- 
sivities are proportional to p 2 , as is the case 
of the recombination lines in plasma. The 
latter regime modifies the analysis. In par- 
ticular, for the shallow spect rum of density, 
Chep urnov & Lazarian I (12006) showed that the 
power spectrum of density P p ~ k~ a has a shal- 
low spectral index a < 3 emissivity spectrum 
P e ~ k" 1 is a e - 2a-3 and this index should be 
used in all the expressions obtained of the VCA 
and VCS techniques. For the steep power law 
index of density, the power law indexes of the 
emissivity and density coinside for sufficiently 
large wavenumbers k. 




Sweet-Parker model 





Turbulent model 



blow up 

Fig. 2. Upper plot: Sweet-Parker model of recon- 
nection. The outflow is limited by a thin slot A, 
which is determined by Ohmic diffusivity. The other 
scale is an astrophysical scale L » A. Middle 
plot: Reconnection of weakly stochastic magnetic 
field according to LV99. The model accounts for 
the stochasticity of magnetic field lines. The out- 
flow is limited by the diffusion of magnetic field 
lines, which depends on field line stochasticity. Low 
plot: An individual small scale reconnection re- 
gion. The reconnection over small patches of mag- 
netic field determines the local reconnection rate. 
The global reconnection rate is substantially larger 
as many independent patches come together. From 
lLazarian et al. I fc004h . 



3. Reconnection and Reconnection 
Diffusion 

It is generally believed that magnetic field em- 
bedded in a highly conductive fluid preserves 
its topology for all tim e due to magne t ic fields 
being frozen-in (e.g., lAlfven 1 1 1942k iParker I 
1979). Although ionized astrophysical objects 
are almost perfectly conducting, they show 
indications of changes in topology, "mag- 
netic reconnection", on dynamical time scale s 
(e.g., lLovelace|[l976l: iPriest & Forbes] l2000l) . 
Reconnection can be observed directly in the 
solar corona (e.g.. | Yokovama & Shi bata 1995; 
iMasuda et al. Ifl994l) . but can also be inferred 
from the existence of large-scale dynamo activ - 
ity inside stellar interiors (e.g.. |Parkerlll993l) . 



Also Solar fl ares (Sturrock 1966|) and y-rav 
bursts (e.g., iFox et al.T 20051 : Galama et al. I 
1998) are usually associated with magnetic re- 
connection. 

To understand the difference between re- 
connection in astrophysical situations and in 
numerical simulations, one should recall that 
the dimensionless combination that controls 
the reconnection rate is the Lundquist num- 
befl defined as S — L x Va/<1, where L x is 



5 The magnetic Reynolds number, which is the 
ratio of the magnetic field decay time to the eddy 
turnover time, is defined using the injection veloc- 
ity V; as a characteristic speed instead of the Alfven 
speed Va, which is taken in the Lundquist number. 
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the length of the reconnection layer, Va is the 
Alfven velocity, and A = t]c 2 /4jt is Ohmic 
diffusivity. Because of the huge astrophysi- 
cal length-scales L x involved, the astrophysi- 
cal Lundquist numbers are also huge, e.g. for 
the IGM they can be as high as 10 20 , while 
present-day MHD simulations correspond to 
S < 10 4 . As the numerical efforts scale as 
h\, where L x is the size of the box, it is feasi- 
ble neither at present nor in the foreseeable fu- 
ture to have simulations with sufficiently high 
Lundquist numbers. 

Observations have always been suggestive 
that magnetic reconnection can happen at high 
speed, in spite of theoretical difficulties to ex- 
plain the effect. At the same, the phenomenon 
of solar flares was suggestive that magnetic 
reconnection may be slow in order to ensure 
the accumulation of magnetic flux and sud- 
denly gets fast to explain the observed fast re- 
lease of energy. A model that can naturally ex- 
plain this and other observational manifesta- 
tio ns of magnetic reconnection was proposed 
in lLazarian & Vishniac I d 19991) (LV99). The 
model appeals to the ubiquitous astrophysical 
turbulence as a universal trigger and controller 
of fast reconnection. 

To deal with strong, dynamically important 
magnetic fields, LV99 proposed a model of fast 
reconnection in the presence of sub-Alfvenic 
turbulence (see Figure They identified 
stochastic wandering of magnetic field-lines as 
the most critical property of MHD turbulence 
which permits fast reconnection. As we dis- 
cuss more below, this line-wandering widens 
the outflow region and alleviates the control- 
ling constraint of mass conservation. The LV99 
m odel has been succ essfully tested recently 
in iKowal et al. I d2009h (see also highe r reso- 
lution results in lLazarian et al. I d2010h ). The 
model is radically different from its predeces- 
sors which also appealed t o the effects of turbu- 
lence. For instance, unlike I Speiser I dl970D and 
Ijacobson & Mosesl d 19841) the model does not 
appeal to changes of microscopic properties of 
plasma^. 



6 The nearest progenit o r to LV99 w as the work 
of iMatthaeus&Lamkinl dl985l Il986h . who stud- 
ied the problem numerically in 2D MHD and 



The LV99 model justifies the notion of tur- 
bulent mixing perpendicular to magnetic field 
lines. Indeed, LV99 showed that the GS95 
model gets self-consistent only in the presence 
of the turbulence-induced reconnection with 
the rates predicted in LV99. Otherwise, the for- 
mation of the magnetic knots would change the 
character of the turbulent interactions. 

The understanding of fast magnetic re- 
connection in the presence of turbulence in- 
duced the notion of "r econnection diff usion" 
that was described in Lazarian (2005|) and 
later used for describing different phenomena 
from star formation to h eating of magnetic 
filaments in IGM (e.g., ISantos-Lima et afTI 
2010; Lazaria n et al. 2010). The same concep t 
was implicitly used earlier in lCho et"al. I d2003l) 
where it was claimed that the heat conductiv- 
ity of the IGM is influenced by the heat ad- 
vection by turb ulent eddies. E xplicit calcula- 
tions done by lLazarian ] d2006h show that the 
heat conduction by turbulent eddies mixing 
magnetic field perpendicular to the local direc- 
tion of magnetic field is the dominant way of 
heat transport in clusters of galaxies. The ef- 
fect of reduced mean free path of thermal elec- 
trons induced by turbulence that we discussed 
above (Sect. 2) increases the relative impor- 
tance of thermal transfer via reconnection dif- 
fusion. Rigorous arguments justifying the con- 
cept of reconnectio n diffusion can be found in 
lEvinketal.ld201ll) . 

4. Cosmic ray acceleration 

Radio observations of galaxy clusters probe 
particle acceleration by shocks and turbulence 
in the IGM (Brunetti 201 1, this conference for 
review on physics of cosmic rays (CR) in the 
IGM). In this Section we briefly discuss the 
importance of turbulence in the acceleration of 
CR and the connected issue of CR acceleration 
induced by magnetic reconnection. 

who suggested that magnetic reconnection may be 
fast due to a number of turbulence effects, e.g. 
multiple X points and turbulent EMF. However, 
iMatthaeus & Lamkml < fl98lll986l) did not observe 
the important role of magnetic field-line wandering, 
and did not obtain a quantitative prediction for the 
reconnection rate, as did LV99. 
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4.1. Acceleration by magnetic 
turbulence 

The interaction of turbulence and cosmic rays 
is a vital component of models of CR prop- 
agation and acceleration. It has been a con- 
cern from the very b e ginning of CR research 



(e^ iGinzbir g 1 1 1 966t I Jokipiil Il966t IWentzell 



1969). To account for the interaction properly, 
one must know both the scaling of turbulence, 
the changes with time of turbulence spectrum 
due to the damping processes (e.g. with CR), 
and the interactions of turbulence with various 
waves produced by CRs. 

Clusters of galaxies present magnetic 
fields of the largest extend and they are also 
considered on the role of the accelerators of 
the ultra high energy CR. The acceleration 
of particles in large (Mpc) regions in galaxy 
clusters is generally believed to happen via 
the second order Fermi process as a result 



teractions (e.g., 


Brunetti & Lazarian 1 120071: 


Petrosian & East 


2008: iBrunetti et all 12008). 



Similarly, acceleration by magnetic turbu- 
lence is a very robust process that is likely 
be important for Solar flares, gamma ray 
bursts and ma ny other astrophysical envi- 



ronments (e .g., Hamilton & Petrosian 1992k 



Miller et al 



" 1996; 



_ Schlick eiser & Dermerl 

2000l:lDermer & Humill200ll) . 

MHD turbulence is the most important 
for the acceleration of particles of largest 
energies and it is vital to use the theoret- 
ically justified and numerically tested rela- 
tions in the studies of particle acceleration. 
From the start of the w ork in this direction 
(e.g., IChandranl l2000) it became clear that 
the earlier models for the acceleration and 
propagation of energetic particles that were 
based on ad hoc representation of turbulence 
are in error of many orders of magnitude as 
far as Alfvenic perturbations are concerned. 
lYan & Lazarian I d2002l 2004) identified com- 
pressible fast modes as the principal agent for 
CR acceleration by MHD turbulence. As the 
aforementioned modes, unlike Alfvenic ones, 
are subject to rather strong damping, the de- 
scription o f the acceleration gets more compli- 
cated. In IBrunetti & Lazarianl d2007l) we de- 



rived a comprehensive picture of compressible 
turbulence in galaxy clusters and studies CR 
acceleration considering all the relevant damp- 
ing processes, with the results providing good 
correspondence with observations. More re- 
cently we extended this formalism to the case 
of the reacceleration of CR and of the sec- 
ondary par ticles generated in the IGM via pp 
collisions (Brunetti & Lazarian 1 1201 
In addition, the accuracy of the particle ac- 
celeration using analytical theory has been 
improved by extending the quasi-linear the- 
ory to the regime of substantial perturbations 
of magnetic field and applied to the case of 
Solar flares dYan & Lazarian Il2008t lYan et al71 
2008). The improved theory has been suc- 
cessfully tested with direct tracing of CR 
trajectories in data cubes obtained with re- 
sults of direct M HD simulations of turbulence 
dBeresnvak et al. 1201 ll) . Future applications of 
these extensions to the case of galaxy clusters 
will be important. 

Compressible turbulence interacts both with 
CR and with thermal particles. This interac- 
tion may also induce magnetic field perturba- 
tions (trough plasma instabilities, e.g. Sect. 2) 
that may further come into play in the particle 
acceleration process. First attempts in this di- 
rection suggest that the fraction of turbulence 
that goes into CR acceleration increases when 
turbulen t-induced instabilities are taken into 
account dBrunetti & Lazarian 1 1201 lah . 

4.2. Shock acceleration and turbulence 

Here we focus on the importance of turbulence 
in shock acceleration mechanisms. Shock ac- 
celeration is thought to be one of the princi- 
pal accepted mechanisms of energetic particle 
acceleration. The shock induces compression 
and particles trapped between magnetic fluctu- 
ations ahead and behind shocks fill the accel- 
eration every time they bounce back and forth 
between converging fluctuations. This is an ef- 
ficient way of accelerating particles which re- 
sults in the energy gain per bouncing to in- 
crease as the first power of the ratio of the par- 
ticle velocity to that of light, i.e. v/c, making 
this process known as the first order Fermi ac- 
celeration. 
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Shock acceleration in galaxy clusters is be- 
lieved to contribute the most of the CR (pro- 
tons), while shock acceleration of CR electrons 
is the most pop ular model to expla in the origin 
of radio relics dEnBlin et all (Q998), Ryu 20 1 1 , 
Briiggen 201 1, this conference for review). 

The necessity of particles to bounce back 
and forth limits the efficiency of the accelera- 
tion of high energy particles through a require- 
ment that the energetic particle should have the 
Larmor radius less than size of the magnetic 
fluctuations that they bounce off. Therefore 
to increase the energy of the accelerated par- 
ticles one should have strong magnetic field 
and strong magnetic fluctuations both in the 
preshock and postshock regions. The situation 
with the postshock region is relatively sim- 
ple. Gas passing through the shocks is known 
to cre ate turbulence (e.g jGiacalone & Jokipii 
12001) . The turbulence is known to increase 
the magnetic field energy, enabling particles to 
scatter efficiently and return to the shock re- 
gion for further acceleration. For the preshock 
region, most work was concentrated on in- 
stabilities that can enhance magnetic field. 
The most comm only discuss ed is the so-called 
Bell instability dBell 1 12004 which is a non- 
resonant current driven instability, that can in- 
cr ease magnetic fie l d in f ront of the shock. 
In Beresn vak et al. I (120091) we proposed that 
a turbulent generation of magnetic field is 
happening in front of the shock, in the re- 
gion which is called precursor. The proper- 
ties of precursor and its formation in front of 
the shock are des cribed in the literature (e.g., 
Mal kov&Drurv 11201 ll) . As the precursor in- 
teracts with the density inhomogeneities pre- 
existing in the medium in front of the shock, 
it gets perturbed, creating vorticity and turbu- 
lence. New studies of t urbulent amplifi cation 
of magnetic field (e.g., ICho et al. 2009) pro- 
vide the rates of magnetic field amplification 
by turbulence. These ra tes were made use of in 
iBeresnvak et al. I (120091) to obtain the values of 
the turbulent magnetic field that is generated 
in front of the shock. The corresponding esti- 
mates show that the preshock magnetic fields 
produced via this process are larger than those 
arising from the Bell instability and that they 
account for cosmic ray acceleration in galactic 




v4 

Fig. 3. CR spiral about a reconnected magnetic field 
line and bounce back at points A and B. The re- 
connected regions move towards each other with the 
reconnection velocity Vr. The advection of cosmic 
rays entrained on magnetic field lines happens at the 
outflow velocity, which is in most cases of the order 
of Va- Bouncing at points A and B happens because 
either of streaming instability induced by energetic 
particles or magnetic turbulence in the reconnection 
region. In reality, the outflow region gets filled in 
by the oppositely moving tubes of reconnected flux 
which collide only to repeat on a smaller scale the 
pattern of the larger scale reconnection. Thus our 
cartoon also illustrates the particle acceleration tak- 
ing place at smaller scales. 

supernovae shock up to 10 15 eV, the so-called 
"knee" of the cosmic ray spectrum. 

Further development of this direction 
presents a very promising avenue of the cos- 
mic ray acceleration research. Interesting boot 
strap processes are likely to be at work as gen- 
eration of magnetic fluctuations in front of the 
shock increases the efficiency of the acceler- 
ation, contributing to the development of the 
precursor. 

4.3. Acceleration induced by magnetic 
reconnection 

An important consequence of fast reconnec- 
tion of turbulent magnetic fields that we 
discussed in Sect. 3 is the formation of a 
thick volume filled with reconnected magnetic 
flux loops. These 3D loops contract, present- 
ing favorable conditions for energetic parti- 
cle acceleration. This process of first order 
Fermi acceleration of energetic particles in 
reconnection regions has been described in 
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Ide Gouveia Dal Pino & Lazarian ] (120051) (see 
also Fig. [3]) for the situation when there is no 
back reaction of the accele rated particles on the 
reconnected magnetic flux. lDrake et al. I(l2006l) 
appealed to a similar process within their pre- 
ferred model of collisionless reconnection and 
proposed that firehose instability can play a 
role of the feedback for the accelerated parti- 
cles. 



More recently, the acceleration in recon- 
nection regions has obta ined observational 
suppor t. It was suggested in Lazarian & Opher 
(2009) that anomalous CR measured by 
Voyagers are, in fact accelerated in the re- 
connection r e gions of magnetopose (see also 
iDrake et al. I ( 20101) '). Such a model explains 
why Voyagers did not see any signatures 
of acceleration passing the Solar system ter- 
mination shock. In a se parate development, 
lLazarian & Desiatil (1201 Oh appealed to the en- 
ergetic particle acceleration in the wake pro- 
duced as the Solar system moves through in- 
terstellar gas to explain the excess of cosmic 
rays of the range of both sub-Tev and multi- 
TeV energies in the direction of the magne- 
totail Magnetic reconnection is ubiquitous in 
astrophysical circumstances and therefore it 
is expected to induce acceleration of parti- 
cles in a wide range of astrophysical environ- 
ments. For instance, the process has been al- 
ready discussed for the accel eration of par- 
ticles in gamma ray burs ts (lLazarian et al. I 
12003b IZhang & Yan 1 120111) and microquasars 
dde Gouveia Dal Pino & Lazarian 1 120051) . We 
expect the process to be important for the ac- 
celeration of protons and electrons in galaxy 
clusters. 



Numerical 2D s imulations presented 
in IDrake et al. I d2010h confirmed high effi- 
ciency of particle acceleration in regions of 
magnetic reconnectio n. However, results in 
Lazarian e t al. I d2010l) show that the process of 
acceleration happens rather differently in 2D 
and 3D situations. The 3D geometry shows a 
wider variety of acceleration regimes and this 
calls for much more detailed studies of the 
acceleration. 



5. Summary 

The main points of our review can be summa- 
rized as follows 

- Turbulence is essential for understanding 
of the IGM. On large scale the descrip- 
tion of turbulence obtained in MHD can be 
used. Compressions induced by turbulence 
induce instabilities in the IGM, changing 
the mean free path of thermal ions. This 
should extend the range over which the 
MHD description of turbulence is applica- 
ble. 

- Studies of turbulence in the IGM can get a 
boost if Doppler-broaderned spectral emis- 
sion and absorption lines are used. The 
techniques originally developed and suc- 
cessfully used in the interstellar research, 
namely Velocity Channel Analysis (VCA) 
and Velocity Correlation Spectrum (VCS) 
are promissing for studing of turbulence in 
the IGM. 

- Magnetic reconnection happens fast in tur- 
bulent media, which makes the models of 
MHD turbulence self-consistent. Fast mag- 
netic reconnection makes MHD turbulence 
somewhat similar to hydrodynamic if one 
considers turbulent motions perpendicular 
to the local direction of magnetic field. 
Such motions can induce a process of "re- 
connection diffusion" which efficient heat 
transfer in the IGM. 

- Magnetic turbulence is very important for 
particle acceleration in clusters of galaxies. 
It can accelerate particles through direct 
interactions with turbulent fluctuations. 
However, it can also modify shocks, in- 
ducing magnetic field generation in shock 
precursors and increasing the efficiency 
of high energy particle acceleration by 
shocks. In addition, it can enable fast mag- 
netic reconnection which can accelerate 
particles within the thick reconnection re- 
gions. 
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